High Performance Computing for Solving Fractional Differential Equations with Applications

Fractional calculus is the generalization of integer-order calculus to rational order. This subject has at least three hundred years of history. However, it was traditionally regarded as a pure mathematical field and lacked real world applications for a very long time. In recent decades, fractional calculus has re-attracted the attention of scientists and engineers. For example, many researchers have found that fractional calculus is a useful tool for describing hereditary materials and processes. It has been used to model the properties of viscoelastic materials and anomalous diffusion. Other applications of fractional calculus include signal processing, control of dynamic system, fractal theory, finance. In this thesis, we have investigated several applications of fractional calculus and the use of multi-core hardware architecture for solving fractional differential equations. Within heat theory, we have studied fractional generalized Cattaneo equations and pointed out that the fractional heat equations may give negative absolute temperatures. Related to elastography, we have investigated the use of a fractional wave equation to describe the shear propagation induced by radiation force. We have concluded that there is a possibility of biased estimation of shear modulus. Numerical simulation of fractional partial differential equations is a time-consuming task due to the non-local property of fractional derivatives. We have shown that optimization techniques and parallel computing can reduce the long simulation time. We have also developed performance models which can give deep understanding of the optimization techniques and predict the simulation time of both serial and parallel implementations. Last but not least, we have demonstrated that parallel solvers of three-dimensional time-fractional diffusion problems are well suited for cutting-edge parallel hardware.